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ABSTRACT 

We present CO observations of nine ULIRGs at z ~ 2 with /jy(24/im) > 1 mJy, previously confirmed 
with the mid-IR spectra in the Spitzer First Look Survey. Ah targets are required to have accurate 
redshifts from Keck/GEMINI near-IR spectra. Using the Plateau de Bure millimeter- wave Interfer- 
ometer (PdBI) at IRAM, we detect CO J(3-2) [7 objects] or J(2-l) [1 object] line emission from eight 
sources with integrated intensities {^-9)a. The CO detected sources have a variety of mid-IR 
spectra, including strong PAH, deep silicate absorption and power-law continuum, implying that these 
molecular gas rich objects at z ~ 2 could be either starbursts or dust obscured AGNs. The measured 
line luminosity "^00(3-2) (1.28-3.77) x 10^^ K km/s pc^. The averaged molecular gas mass M^^ is 

1.7xlO^°M0, assuming CO-to-H2 conversion factor of 0.8 Mo/pkm/spc^]. Three sources (33%) - 
MIPS506, MIPS16144 & MIPS8342 - have double peak velocity profiles. The CO double peaks in 
MIPS506 and MIPS16144 show spatial separations of 45kpc and 10.9 kpc, allowing the estimates of 
the dynamical masses of 3.2x10^^ sin~^(i) M© and 5.4x 10^^ sin~^(i) M© respectively. The implied 
gas fraction, Mgas/Mdyn, is 3% and 4%, assuming an average inclination angle. Finally, the analysis 
of the HST/NICMOS images, mid-IR spectra and IR SED revealed that most of our sources are 
mergers, containing dust obscured AGNs dominating the luminosities at (3-6) /im. Together, these 
results provide some evidence suggesting SMGs, bright 24 /im ULIRGs and QSOs could represent 
three different stages of a single evolutionary sequence, however, a complete physical model would 
require much more data, especially high spatial resolution spectroscopy. 

Subject headings: galaxies: infrared luminous - galaxies: starburst - galaxies: high-redshifts - galaxies: 
evolution 



1. INTRODUCTION 

Molecular gas holds one of the keys to fundamental 
questions about the formation and evolution of galax- 
ies. Particularly, at high redshift, significant galaxy 
growth in its stellar population and black hole occurs 
during dusty, gas rich phases. The importance of infrared 
(IR) luminous, gas rich populations has been high-lighted 
by observations from the Infrared Space Observatory 
(750), ground-based (sub) millimeter cameras, and more 
recently, the Spitzer Sp ace Telescope ([Genzel fc Cesarsk^ 



2QQQl|Blain et al.||2QQ2| |Lagache et al.||2UU5| l^oiter et a. 
2QQ8p . Particularly, Spitzer dramatically improved our 
ability to probe high redshift (z > 1), dusty galaxies down 
to Ljji^IO^^Lq^ increasing the number of detections 
of such sources by orders of magnitude in comparison 
with ~ (200-300) z ^2 gas rich, sub- millimeter galaxies. 
Several recent studies using the Spitzer InfraRed Spec- 
trograph (IRS; jHouck et al.|[2004) have confirmed with 
mid-IR spectra a population of 24 /im-selected, Ultra- 
luminous infrared galaxies (ULIRGs) at z ^ 2 (Yan et al. 
2005] IHouck et al.||2QQ5l [Weedman et al.||2UU6[ |Yan et 



al.| l2007[ ISajina et al.||2008| |Dey et al.||20081 IDasyra 
al. |2009t. One important feature of these studies is tf 



et 
that 



the samples are selected to have 24 /im fluxes brighter 
than ~ ImJy, placing them among the most I R lumi- 
nous objects dit z ^2 with Ljr^ 10^'^-^~^^-^Lq (Sajina 
et al.|[2QQ7 , 2008). Other gas rich galaxies with compa- 



rable or more Ljr at z > 2 are sub- millimeter galaxies 
(SMGs) and optically bright QSOs with CO detections. 
By definition, SMGs are selected by their far-infrared 
emission, produced by larger, colder dust grains, in con- 
trast, Spitzer bright 24 /im ULIRGs are selected in the 
mid-IR, probing smaller and hotter dust grains. 

SMGs, Spitzer bright 24 /im ULIRGs and gas rich 
QSOs are all high-z systems with extremely high bolo- 
metric luminosities, leading to the important question 
of how these systems are related or different. Studies 
of the gas content, particularly CO gas , of SMGs have 



made significant progress in recent years (Neri et al. 2003 
Greve et al. 2005; Tacconi et al. 2006 2008). The gen- 
eral consensus for this population is that SMGs contain 
a reservoir of 10^^~^^ Mq of molecular gas, distributed 
over a small area of i?i/2^2kpc. This large amount 
of molecular gas is being converted to stars at a rate 
of ~ lOOOM0/yr, typically over several dynamical time 
scales of ~ lO^yrs. Recent several studies have also ex- 
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amined the mid-IR spectral properties of SMGs (|Lutz et 
al.|2005|[Valiante et al.|2'QQ7l|PQpe et al.|2008||Menendez- 
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TABLE 1 
Sample Targets 



ID 


RA 

deg 


DEC 
deg 


z(Opt)'' 


z(CO) 


MIPS429 


259.04930 


59.20366 


2.213 


2.2010 


MIPS506 


257.91080 


58.64405 


2.470 


2.4704 


MIPS8196 


258.79282 


60.16533 


2.586 


b 


MIPS8327 


258.89908 


60.47375 


2.441 


2.4421 


MIPS8342 


258.54813 


60.18591 


1.562 


1.5619 


MIPS15949 


260.28842 


60.25035 


2.122 


2.1194 


MIPS16059 


261.11850 


60.25922 


2.326 


2.3256 


MIPS16080 


259.68655 


60.02107 


2.007 


2.0063 


MIPS16144 


261.09207 


59.53077 


2.131 


2.1280 



^ The optic al spe ctroscopic redshifts are taken from 
[Sajina et al.) ( |2008| >. 

" MIPS8196 has no significant CO hne detection, thus, 
its CO redshift is not determined. 



Delmestre et al.||2QQ7l |2QQ9[ ). The one clear difference is 
that the observed 24 /im tluxes of SMG are on average 
several 100 /iJy, or a factor of (2-3) smaller than their 
Spitzer -selected counterparts. This is consistent with 
the conclusion that bright 24 /im z ^2 ULIRGs are star- 
burst/ AGN composite systems, with more AGN heat 
ing the hot dust grains, thus generating elevated contin 
uum emissi on at (3 - 6) iim than that of SMGs (i Sajina 
et al.||2008| iPolletta et al.||2008|). In contrast, U 



Polletta et al.||2008| ). In contrast, tor the 
majority of'BMGs, star burst dominates the IR lumi- 
nosity, while AGN contributes only a small fr action of 
LiR ( Valiante et aO2007 [[Pope et al .||2008| [Alexander et 
[al. 2008; Menendez-D elmestre eT^ . 2009). In addition, 
abundant CO molecular gas has been detected among 
optically selected, high-z QSOs and ra dio galaxies, espe- 
cially gravitationally lensed systems. S olomon fc Vander? 
[Bout (2005) reviewed and compiled a list of 23 QSOs and 
radio galaxies whi ch have at least low resolution CO line 
observations (also |Gr eve et al.' (2005) for a similar list). 

In this paper, we will present the measurement of cold 
molecular gas masses among bright 24 jam ULIRGs at 
z ~ 2, determine their gas dynamics, and understand how 
they differ from SMGs and gas rich QSOs & radio galax- 
ies. To directly address these questions, we obtained CO 
interferometric observations for a sample of bright 24 /im 
ULIRG s at 2; ^ 2 selected from our previo usly published 
studies ( Yan et al.|2007 Sajina et al.|2008' ). This paper is 
organized as follows. § 2 describes the sample selection, 
ancillary data used for the analyses in this paper, and 
the CO interferometric observations and data reduction. 
§ 3 presents the results, and discusses the physical impli- 
cations by combining all available imaging and spectro- 
scopic data. To understand how 24 jam ULIRGs differs 
from SMGs and high-z QSOs, we make comparisons to 



the QSOs and radio galaxy sample compiled by [Solomon 
& Vanden Bout ( 2005 ) and also to a small subset of QSOs 



which have reliable dynamic masses from high resolution 
CO observations and robust black hole mass estimates 
from UV spectroscopy. § 4 summarizes the main conclu- 
sions from this paper. Throughout the paper, we adopt 
an l^M =0.27, I^a =0.73 an d Hq = 71 kms"^ Mpc"^ cos- 



mology dSpergel et aL||2003 ) 



2. OBSERVATIONAL DATA 



The ground transition temperature for CO J (1-0) is 
only 5.5 K, in contrast for H2 500 K. Therefore, CO ro- 



tational emission lines are easily produced by collisional 
excitation. It has been shown that CO flux Ico linearly 
correlates with the column density of molecular hydro- 
gen H2, providing a sensitive t racer of bulk of the cold 
molecular gas in the Unive rse (Young & Scoville 1982^ 
[PTckman et al. 1986 Solomon et al. 1987). Besides the 
molecular gas content of a galaxy, CO interferometric ob- 
servations can also probe the spatial distribution as well 
as velocity field of molecular gas. Therefore, we obtained 
new CO emission line observations for nine bright 24 /im 
ULIRGs at z - 2 using the PdBI at IRAM. In this sec- 
tion, we describe the selection of our CO targets, and 
summarize the CO observations and all of the ancillary 
all of the data. 

2.1. The sample and ancillary data 

The nine targets observed by the PdBI were selected 
from a large sample of z ^ 1-2 ULIRGs w ith Spitzer 
mid-IR spectra published in | Yan et al. (2007). This par- 
ent sample consists of 52 sources mitially selected with 
24 /im flux density brighter than 0.9 mJy, and very red 
24-to-8 /im and 24-to-R colors from the 4 square degree 
Spitzer Extragalactic First Look Survey (XFLSj^ The 
subsequently obtained Spitzer mid-IR spectra covering 
the rest-frame 4-20 /im determined that 74% of the sam- 
ple is at 1.5<z< 3.2, confirming the effectiveness of the 
initial color selection (Yan et al. 2004| ). The total infrared 
luminosities (L3_iooo/xm) are in the range of lO^^^^^L©. 
The complete mid-IR spectral and SED analysis shows 
that at least '^75% of the sample contain dust obscured 
AGNs, which dominate the mid-IR (3-6) jam luminosi- 
ties but star for mation still contribute mo st of the far- 
IR luminosities ([ Sajina et a'L] |2007[ |2008[). Of the 52 
sources, 44 galaxies were observed at 1.2 mm with the 
IRAM 30 meter telescop e using the 11 7 element version 
of the MAMBO array ( [Kreysa et aL| ri998). Of these, 
7 are detected at > 3cr (14 at 2cr), with an average 
rms ~ 0.6 mJy. T he 1.2 mm data of t he fu h sample has 
been published in Lutz et al. (2005) and Sajina et al. 
(2008). 

To measure cold molecular gas content and to quantify 
how it changes with mid-IR properties among these z^2 
Spitzer ULIRGs, we obtained CO interferometric obser- 
vations for sources with a broad range of mid-IR spectral 
properties, including sources with strong PAH emission, 
deep silicate absorption, and mid-IR power-law continua. 
The primary selection of our CO targets was the avail- 
ability of accurate spectroscopic redshifts from Keck and 
GEMINI. This criterion has limited our targets to a small 
subset of the full sample. All of our targets have 1.2 mm 
observations from MAMBO. Of the nine sources, only 
one have 1.2 mm SNR > 3cr, seven with SNR -1.81 - 2.58, 
and one with SNR< 1.5. Excluding the > 3a source, we 
stacked the remaining 8 sources, yielding an averaged 
1.2 mm flux of 1.06 ± 0.18 mJy. This value is higher than 
the averaged value of 0.5±0.1mJy for the 44 sources 
observed at 1.2 mm in the original Spitzer sample (52 
objects) (Sajina et al. 2008). This implies that our CO 
targets systematically have more far-IR emission, in com- 
parison with the full Spitzer mid-IR spectroscopic sample 
at z^2. 

1 i?(24, 8) = logio(^/^(24/im)/z/^(8/im) > 0.5; 
R{24, 0.7) = logio(^/424Mm)/^/.(0.7/im) > 1.0 



Detections of CO Molecular Gas in Spitzer z ~ 2 ULIRGs 3 



TABLE 2 

CO OBSERVATIONS 



ID 


Lines 


^obs 


Config. 


Beam^ 


PA Time^ 


Noise/chan. 


Chan. width 


RMS^ 






GHz 






hours 


mJy/chan 


km/s 


km/s mJy/beam 



MIPS429 


(3-2) 


107.642 


C 


4.8" X 3.89" 


68° 


6 


0.4 


83 


664 


0.2 


MIPS506 


(3-2) 


99.679 


D 


5.02" X 4.37" 


71° 


4 


0.35 


75 


525 


0.25 


MIPS8196 


(3-2) 


96.456 


C 


5.6" X 4.18" 


65° 


4.5 


0.39 


62 






MIPS8327 


(3-2) 


100.464 


C 


2.86" X 2.5" 


71° 


9 


0.7 


60 


540 


0.12 


MIPS8342 


(2-1) 


89.913 


D 


6.35" X 4.15" 


56° 


9.5 


0.6 


67 


670 


0.1 


MIPS15949 


(3-2) 


111.046 


C 


2.63" X 2.3" 


69° 


10.1 


0.4 


54 


756 


0.12 


MIPS16059 


(3-2) 


103.999 


C 


4.77" X 3.83" 


89° 


5.5 


0.5 


87 


522 


0.25 


MIPS16080 


(3-2) 


114.997 


C 


4.61" X 3.52" 


75° 


10 


0.7 


72 


360 


0.35 


MIPS16144 


(3-2) 


110.407 


C 


2.44" X 2.15" 


62° 


5.4 


0.5 


81 


810 


0.13 



^ Beam is defined as the semi-major axis times the semi-minor axis in arcseconds. 
^ This is the total on-target integration time, calculated equivalent to 6 antenna. 

^ The summed CO maps are integrated over a velocity width A with the center on the peak of the CO emission line. 
^ This noise is RMS for the summed CO map made with the listed parameters. 



TableJTjand Table|2]list the coordinates, redshifts (from 
both C(J and near-lK spectra) and the CO observational 
parameters for the 9 sources observed by PdBI. Table [3] 
shows the ancillary observations from optical, near- IR, 



Sajina et aL[2QQ8[ 
3 obtained usinff 



for 



far-IR and 20 cm photometry (see 
details). The near-IR spectra were 
SPEC at the Keck and NIRI at the GEMINI. We have 
high resolution H-band images from HST mCMO S. This 
dataset has been published in Dasyra et al. (2008) and is 
included in the analysis of a larger sample ot z ^0.3-2 
bright 24 /im ULIRGs in (Zamojsk i et al.,t2QQ9l ). 

2.2. CO emission line observations 

The observations were carried out in 2007-2008 with 
the PdBI in D (2 objects) and C (7 objects) configu- 
ration in order to maximize the sensitivity. The 3-mm 
receivers were tuned to the central observed frequency 
according to the Near-IR spectroscopic redshifts. Of the 
total nine sources, we used the PdBI to observe the CO 
J(3-2) transition (z/^est = 345.796 GHz) for eight objects 
and the J(2-l) transition (z/^est = 230.538 GHz) for one 
object at redshift of 1.562. Typically, our sources were 
observed under good weather conditions with 1-2 tracks 
for a total of (5 -10) hours, 6 antenna equivalent, on- 
target integration. Table [2] lists the observation parame- 
ters. The data reduction takes places in two stages using 
the IRAM GILD AS software. First, the raw data are cal- 
ibrated using the software CLIC, developed at IRAlVj^ 
Care is given to the phase, fiux and amplitude calibra- 
tions by rejecting any calibration anomalies due to me- 
teorological conditions or electronics. After calibration, 
CLIC generates the visibility data, ie. "uv" tables. Sec- 
ond, based on these uv tables, we extract CO maps and 
spectra using the IRAM MAPPING softwar^ The syn- 
thesized, clean beam size is typically elongated, roughly 
(2.4^x2.4") to (5.6^x4.18") for the C configuration and 
(5.02"x4.37") to (6.2"x4.15") for the D configuration. 



complete documentation 



available 



| http://ww w.iram.fr/IRAMFR/GILDAS/doc/html /clic- 
Lhtml/clichtml 

I'he documenta tion is at 

■htt p : / / www . ira m.fr/IRAMFR/GILDAS/doc/html/map- 
html/ map. html 



at 



] 



The corresponding spatial resolution in linear size ranges 
from 20.5 to 53kpc. For comparison, for the Spitzer 
24 /im band, the Fuh Width Half Maximum (FWHM) 

of a point source is 6 and the NICMOS H-band PSF is 
0.22". 

Figure [1] shows the integrated CO maps, with the small 
insert at the left bottom of each panel indicating the 
beam size and positional angle. Before producing the 
summed CO maps, we rebinned the data with an orig- 
inal 10 MHz channel width by a factor of (3-4), yield- 
ing smoothed spectral data cubes with a channel width 
roughly (54-83) km/s. The CO maps shown in Fig- 
ure [1] were made by summing over a number of veloc- 
ity channels, ranging from 5-14 (equivalently, velocity 
width A ~ 360-810km/s in Table |2|, centered on the 
peak CO line emission. The final maps have la noise 
values ranging from 0.1 to 0.35 mJy/beam. Table [2] lists 
the noise parameters and the parameters used for mak- 
ing the summed CO maps. In Figure [ij the first CO 
contour starts roughly at la with a step size of la for all 
panels. The summed CO line emissions are detected at 
(4.7-9.7)cr for 8 sources. Table[4]lists the integrated CO 
line fiuxes, velocity width over which the CO line inte- 
gration is done, and derived CO luminosity and H2(+He) 
molecular masses. 

The CO line luminosity can be derived from the line 
intensity, following Solomon et al. (1997): 



Leo [i^o] = 1.04 X 10- 



and 



Ico 



^obs 



Dt 



Jykms-iy VGHzV VMpc) 



(1) 



L^o [K km/s pc'] =3.25 x 10%o^obs"'DL'(l + z)-^ 

(2) 

Here L^q is the CO line luminosity in units of 
Kkm/spc^, Ico is the integrated line intensity in 
Jykm/s, Dl is the luminosity distance in units of Mpc, 
Uobs is the observed frequency in GHz. The H2 (+He) 
molecular gas mass can be estimated using the following 
equation: 
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5 10 0.001 

Jy km/s 



Fig. 1. — The integrated CO maps for the 8 detected sources. The smah insert in each panel shows the size and shape of the beam. The 
cross marks the center of 24/im source. The first contour starts at l.lTcr, l.lScr, l.OScr, l.OScr, 1.04cr, l.lScr, Icr, and l.OScr respectively 
for MIPS429, MIPS506, MIPS8327, MIPS8342, MIPS15949, MIPS16059, MIPS16080 and MIPS16144. The contour step size is la for ah 
maps. The rms a (mJy/beam), the channel width and the number of channels over which the CO maps are integrated, are given in Table [2] 



Mr. 



^CO(3-2)/^CO(l-0) 



1 



(3) 



Here a is the CO J (1-0) luminosity to molecular gas 
mass conversion factor in units of [K km/s pc^]~^ and 
Too is the effective brightness temperature of a CO tran- 



sition. Here we assume T( 



CO(3-2) 



/^COj 



1-0) 



1 and 



a = 0.8 McT) (K km /s pc^)~^ ( [Solomon et al.|l997||Downes 
|& Solomon'"1998), the same assumption used tor SlVlGs 
and gSOs (Greve et al.. 2005, Tacconi et al.^ 2006, 2008} . 



Although a more realistic T co(3-2)/Tco(i-o) rat io is 
probably between (0.5-0.7) ( [Wilson et aT] |2008D , we 
adopt our assumption in order to make direct compari- 
son without additional corrections to the published SMC 
and QSO data. 

2.3. Modeling the IR SED 

To derive the total and far-IR luminosities 
(LiR = L3_iooo/im,LFiR = L40-i20Mm), wc fit the spectral 
energy distribution of our sources from near-IR to far-IR 
(see Table [3| . The longest wavelength data available is 
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MAMBO 1.2 mm, which is crucial for constraining the 
far-IR emission. All of the nine targets except one have 
weak 1.2 mm fluxes at 3a or less. For these 8 sources, 
the sigma weighted, averaged flux is 1.06 ± 0.18 mJy 
(5.9<j), providing constraining power in the SED fitting 
in Figure [2] (cyan, open stars). Our SED fitting uses 
three components: pu re star burst template (th e aver- 
aged SMG SED from [Pope et al.| ( [20061 [2008^ ^ 
dot- dashed l ine), reddened type-1 quasar •EED lRichard^ 



et al. (2006), blue dashed line), and a 2Gyr old S ingle 



Stellar Population (SSP) model from Maraston (2005) to 
account for the near-IR emission (purple dot-dot-dashed 
line). The solid red line is the sum of these three 
components. The type I quasar SED (blue) has to be 
reddened substantially (Ay > 5) assuming a MW-type 
extinction (Draine 2003). For most sources this leads 
to a reasonably good fit to the mid-IR continuum 
and silicate absorption feature. However, two of our 
sources (MIPS8342 and MIPS15949) show red mid-IR 
continua but no significant silicate absorption. The 
above approach cannot fit this type of spectra. For these 
sources, we revert to our power-law continuum approach 
as used i n ^Sajina et al. (2008 ). 

At the far- IK, lour sources are detected in the 70 jam 
band, and only one at 160 jam. With so few far-IR 
detections, is the starburst-component really required? 
The only sources for which the answer is obviously 
yes are MIPS8342 (based on its 160 /im detection) and 
MIPS 16144 (as it is a strong- PAH/mm-bright source). 
For MIPS 16080, its 70 /im detection is difficult to explain 
without such a component, unless we assume a somewhat 
cooler AGN template. In particular, the stacked 1.2 mm 
detection (1.06 ± 0.18 mJy) requires a minimum amount 
of star formation, as indicated by the quiescent spiral 
galaxy template shown as the solid pink line in Figure [2] 
(this is a Sd template from Polletta et al. (2007)). 

This naturally begs the question whether the starburst 
or the quiescent Sd template is best suited to fit the SEDs 
of our sources. The choice will greatly impact the cal- 
culated far-IR luminosity. We argue that in order to 
fit simultaneously the observed 70, 160 /im and 1.2 mm 
fluxes, the starburst template is required. Sd templates 
can explain the stacked flux at 1.2 mm, but fail signifi- 
cantly with 70 /im detections and mid-IR spectra. The 
stacked limit at 70 /im is (2 ± 0.68) mJy, high enough that 
Sd template would fail to fit this limit by a large mar- 
gin. Definitive solution to this problem will soon to be 
available from Herschel data with better spatial resolu- 
tion and sensitivity. The luminosities derived from our 
SED fitting are tabulated in Table |5j We note that Lpm 
is dominated by starburst component for all sources in 
this study. 

3. RESULTS 

3.1. Detections of CO Line Emission 

Of the nine sources observed, eight yielded signif- 
icant detections in the CO J (3-2) or J (2-1) transi- 
tions. The CO J(3-2) spectrum of MIPS16144 also 
shows a weak continuum at the level of 0.5±0.1mJy 
(Aobs = 2.717mm). This observed continuum flux is con- 
sistent with the full-IR SED (see the open triangle sym- 
bol in Figure [2|. Figure [3] presents the CO maps, 
HST/NICMOS Il-band images and mid-IR spectra, one 



row per source. The first column shows the CO contours 
overlaid on the 24 /im image (red, FWHM~ 6 ) and the 

HST NICMOS H-band image (dark, FWHM-0.22"). 
All CO contours start at la with an increment contour 
step size of la. MIPS429 has no NICMOS H-band data, 
and we used the available WFPC2 F814W image instead. 
It has bright 24 /im flux but no optical counterpart in 
both ground-based R image and the HST/WFPC2 im- 
age. The second column shows the HST/NICMOS H- 
band images of our sources, providing morphologies at 
the rest-frame ~ 5500A. The third and fourth columns 
are the CO and mid-IR spectra. The mid-IR spectra for 
MIPS16059 and MIPS16080 have higher SNR, and are 
from an ultra-deep Spitzer observation aiming to detect 
water ice and hydrocarb ons among a samp le of deeply 
embedded z-2 ULIRGs ( |Sajina et aL]|2009D . 

The strength of PAH emission and the 3 - 6 jam hot 
dust continua are indicators of star formation and AGN 
respectively. One interesting question is whether the 
mid-IR star formation indicator is correlated with cold 
dust and molecular gas. Figure |3] shows one obvious re- 
sult, that is the detect ability of CO gas is not strongly 
correlated with the mid-IR spectral properties, with the 
caveat that this finding is based on a small number of 
sources. Of the nine sources, seven have mid-IR spectra 
whose prominent features is the broad emission bump at 
(7-8) /im and/or fully or partially observed silicate ab- 
sorption trough. This type of spectra is very characteris- 
tic for highly obscured 24 jam sources at z ^2, as shown 
by several recent Spitz er observations ( Houck et al.|2005" 



|Yan et al.|2005l |2007|) . Local examples with this type a 



spectrum are IRAS F00183-7111 and NGC4418 ( [Spoon 
,et al.^2004a^b^ ,2001 ) . One characteristics of these deeply 
embedded local ULIRGs is their multi-temperature ISM, 
ranging from hot gas near the cent ral nuclei to cold gas to 
outer region ( [Spoon et al.||2004b ). Our CO observations 
show that six of these seven sources with deep silicate ab- 
sorption have cold molecular gas, one (MIPS8196) does 
not have CO detection at Ico ^ 3cr ~ 0.5 Jy km/s. 

The CO detection rate of our sample is high, likely 
due to the improved sensitivity and larger bandwidth of 
PdBI in recent years. Physically, such high detection rate 
implies that among ImJy or brighter 24 /im ULIRGs at 
z~2, a good fraction of them, including SB/ AGN com- 
posites, deep silicate absorption systems, and AGN dom- 
inated galaxies, could have abundant cold molecular gas. 
For these 24 jam bright ULIRGs with substantial black 
holes in their centers, one important yet unsolved ques- 
tion is how cold molecular gas is exactly turned into stars 
and feeds the growth of black holes. The answers will 
come from future high resolution studies with ALMA. 

3.2. Cold molecular gas fraction and far-IR emission 

Table [4] lists the integrated CO line intensity Ico 
in Jykm/s and luminosity L^q for J(3-2) or J(2-l) 
transition. We emphasize again that to derive cold 
molecular gas mass, we assume T'co(3_2)/T'co(i_o) = 1, 
^co(3-2)/^co(i-o) = 1 ciiid CO J(l-O) luminosity to molec- 
ular gas mass conversion factor a = 0.8 Mq/[K km/s pc^]. 
The eight sources have Mgas (1.02-3.02) xlO^^M©, 
with a median value of 1.7xlO^^M0. In order to com- 
pare our sample with other galaxy populations at z > 2 
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TABLE 3 

Ancillary data and derived parameters for the targets 



ID R H Flux Flux" Flux" Flux Flux Flux 

vega AB 24)Ltm TOyLtm IQO/^m 1.2mm 1.4GHz 650MHz 
mag mag mjy mjy mjy mjy mjy mjy 



MIPS429 


<25.5 




1, 


.10 


<4.5 


<30 


1.03±0.57 


<0.08 




MIPS506 


23.44 


21.66 


1, 


.06 


<4.5 


<30 


1.37±0.53 


0.14 ±0.03 




MIPS8196 


22.45 


18.68 


1, 


.50 


5.0 ± 1.4 


<30 


0.99±0.43 


<0.08 




MIPS8327 


23.27 


21.06 


1. 


.16 


<5.8 


<30 


1.03±0.59 


1.4 ±0.06 


3.44 ±0.14 


MIPS8342 


23.51 


21.54 


1, 


.17 


10.7 ± 1.4 


29±11^ 


0.98±0.52 


0.18 ±0.03 


0.86 ±0.11 


MIPS15949 


22.90 


20.73 


1, 


.50 


7.3±1.6 


<30 


1.24ib0.51 


0.16 ±0.02 




MIPS16059 


23.48 


20.49 


1, 


.29 


< 7.2 


<30 


1.20ib0.66 


0.57 ±0.03 




MIPS16080 


22.86 


20.31 


1, 


.10 


5.2 di 1.7 


<30 


0.69ib0.54 


0.34 ± 0.03 




MIPS16144 


23.10 


20.85 


1, 


.12 


< 3.6 


<30 


2.93ib0.59 


0.12 ±0.03 





^ For none-detections at 70 /im and 160 /im, the fluxes are listed as 3cr. 

^ MIPS8342 70 /im position corresponds to multiple 24 /im sources. Here we deblended the 70 /im 
blob, and list only the partial flux belonging to MIPS8342 with the flux error from the total. See 
Sajina et al. (2008) for detail. 



TABLE 4 
Derived CO parameters 



ID" 


FWHMl^ 


FWHM2^ 




Leo 


LljlO^O)^ 


Mass(H2)^ 




km/s 


km/s 


Jy km/s 


lO^Lo 


K km/s pc"^ 


lO^^M© 


MIPS429 


565±80 


777 


1.00 ±0.17 


3.56 


2.9 


2.32 


MIPS506 


300±80 


330 


0.52 ±0.11 


2.16 


1.63 


1.30 


MIPS506'^^^ 


139±50,88±40 


150,87 


0.52 ±0.11 


2.16 


1.63 


1.30 


MIPS8196 






<0.5 


<2.2 


<1.67 


<1.34 


MIPS8327 


253±50 


303 


0.42 ± 0.07 


1.69 


1.28 


1.02 


MIPS8342 


325±70 


406 


0.65 ± 0.07 


0.88 


2.24 


1.79 


MIPS8342^^^ 


134±40,144±40 


128,130 


0.65 ± 0.07 


0.88 


2.24 


1.79 


MIPS15949 


500±117 


129 


1.07 ±0.11 


3.59 


2.72 


2.18 


MIPS16059 


471±54 


426 


0.56 ±0.11 


2.33 


1.76 


1.41 


MIPS16080 


130±60 


159 


0.85 ±0.15 


2.69 


2.03 


1.62 


MIPS16144 


744±217 


824 


1.45 ±0.27 


4.99 


3.77 


3.02 


MIPS16144^^^ 


404±70,404±70 


405,334 


1.45 ±0.27 


4.99 


3.77 


3.02 



^ The sources with double peak CO spectra have two entries per object in this table. The entry 
marked with dbl has the velocity FWHM fltted with double gaussian proflles, thus the second 
column FWHMl has two values per row. 

^ FWHMl is from the fltting with one or two gaussians, see above note. 

^ FWHM2 is directly measured from the data as the full width at the half peak intensity, without 
assuming the gaussian profile. For the double peak CO spectra, FWHM2 are measured for both 
peaks separately. 

^ Here most CO line intensity and luminosity are for J (3-2) transition, except MIPS8342, which 
is J(2-l). To derive H2 mass, we assume T^o(3-2) /'^co(i-O) — 1 ^l^o Tco(2-i)/'^co(i-0) — 1- 
a is CO J (1-0) luminosity to H2 (±He) molecular mass conversion factor. We adopt 
q; = 0.8 M0/[K km/s pc^] for this paper, which is suitable for SMGs and local ULIRGs. 



TABLE 5 

Derived parameters from the ancillary data for the targets 



ID Logio[Lfj^] Logio[L|,jj^] Re// EW(7.7PAH)r^est rg.s^im LogiofLj'^QjjJ 

L0 L0 kpc /im WHz~^ 



MIPS429 


12.72 


12.50±0.15 




0.41 ±0.1 


>7.3 


< 24.4 


MIPS506 


12.93 


12.70±0.09 


1.97 


0.87 ±0.36 


>6.7 


24.67 ±0.14 


MIPS8196 


12.99 


12.45±0.15 


3.31 


0.31 ±0.08 


1.3 ±0.4 


< 24.5 


MIPS8327 


12.83 


12.47±0.21 


2.07 


< 0.3 


2.4 ±0.7 


25.88 ±0.02 


MIPS8342 


12.55 


12.13±0.30 


2.28 


0.58 ±0.14 


0.2 ±0.4 


24.82 ±0.08 


MIPS15949 


12.90 


12.42±0.18 


2.36 


0.31 ±0.06 


0.0 ±0.1 


24.60 ±0.20 


MIPS16059 


12.87 


12.57±0.15 


2.82 


0.29 ±0.07 


2.7 ±0.8 


25.23 ±0.04 


MIPS16080 


12.71 


12.38±0.15 


2.43 


0.11 ±0.05 


2.1 ±0.5 


24.88 ± 0.05 


MIPS16144 


12.72 


12.61±0.12 


3.24 


2.50 ±0.26 


2.0 ± 1.4 


24.46 ± 0.43 



^ LiK is integrated over 8-1000/im, LpiR is over 40-120/im. 

^ Li.4GHz is the monochromatic luminosity at the rest- frame 1.4GHz, derived based on the observed 
radio fiuxes at 1.4 GHz and 650 MHz. 
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Fig. 2. — Near-to-far-IR spectral energy distributions for the 9 targets. The figure legend shows the starburst, AGN and stellar components. 
The large cyan stars at the observed 1.2 mm mark the averaged fiux value of 1.06 zb 0.18 mJy. The down-ward arrows mark the limits of 
undetected sources. For MIPS16144, the open black triangle indicates the rest-frame continuum fiux observed in the CO J (3-2) spectrum. 
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Fig. 3. — The first column shows the the 24 /im image in pink, HST/NICMOS H-band image in black in its original spatial resolution, 
and CO image in contours. In all panels, the CO contours start at la with a step size of la. For a values, see the text for detail. The 
image size is 12 xl2 . The second column shows the H-band morphology. The third and fourth column show the CO and mid-IR spectra. 
The MIPS 16144 CO spectrum shows a detection of continuum, marked by the horizontal dashed line. 



with CO detection, we compiled a list of SMGs based on the published data (Greve et al. 2QQ5[ [Tacconi et 
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al. 2006 '2008), shown in Table [6j This sample of 14 addition, a sample of 19 QSOs and radio galaxi es with 
ISlVlGs includes only sources with significant detections, low re solution CO observations (Solomon & Vanden Bout 
and their median gas mass Mgas is S.OSxlO^^M©. In 2005) is also used as a QSO comparison sample. As de- 
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Fig. 3. — Of 9 objects observed, this is the only source with no significant detection. 



scribed in § |3.3[ we have also constructed the second QSO 
comparison sample, much smaller subset but with reli- 
able dynamical masses from the recent high resolution 
CO observations. Figure Ul left panel shows L^jq ver- 
sus redshift for the 24 jim. U LIRGs and the comparison 
samples of SMGs, QSOs and radio galaxies. The right 
panel in Figure Ul presents the histogram of L^q with 
the dashed lines marking the median values of the three 
populations. The average CO mass of our 24 /im ULIRG 
sample is about a factor of 2 less than the average val- 
ues calculated from the specific comparison samples of 
SMGs and gas rich QSOs currently available in the lit- 
erature. Kolmogorov-Smirnov tests on 24 /im ULIRGs 
versus SMGs and 24 /im ULIRGs versus QSOs derived 
the probabilities of two datasets drawn from the same 
parent sample of 0.19 and 0.189 respectively. We em- 
phasize that the statistics is still too small to draw any 
definitive conclusions for the overall populations of these 
three types of sources. 
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Fig. 4. — This compares the CO integrated hne luminosity as a 
function of redshift for various sources at high- 2;. The SMG data 
are from Greve et al. (2005) and Tacconi et al. (2006). The QSO 

radio galaxy data are from the compilation of Solomon & vanden 

Bout (2005). The dashed lines in the L^q distribution plot (right 
panel) mark the median values for the three galaxy populations. 
The comparison reveals tentative evidence that on average 24 /im 
ULIRGs may have less cold molecular gas than those of SMGs and 
QSOs. 

If far-IR emission is considered mostly from massive 



young stars born in recent starbursts, LpiR can be di- 
rectly related to star formation rate (SFR). Therefore, 

the ratio of Lfir/L^q can be interpreted as describe 
(1) gas depletion time scale, tsf = Mgas/SFR or (2) star 
formation efficiency, SFE = Lpm/Mgas, ^.e. how much 
luminosity can be produced by 1 Mq of gas. If we 
take the SFR[M0/yr] = c * {^^mt:)^ where c = 0.8-3. 



Meurer et al.|T997| . We adopt c = l.S because we use 



depending on the definition of LpiR jKennicutt] [1998 



I^FIR = I^40-i20/xm and the commonly adopted relation 
by Kennicutt (1998) has c = 1.72 with LpiR = Lg-iooo/im- 
The Lftp is measured from the full SED fitting described 



m 



2.3 



The inferred SFR are 32- 787 M© yr-^ with 

the mean value of 585 Mq yr~^. These values are sub- 
stantially larger than normal star forming galaxies, but 
a factor of 2-3 less than what has been observed for 
SMGs at the similar redshifts. Furthermore, if we as- 
sume that 70% of LpiR from SMGs is due t o star for- 
matio n, as su^^ested by some recent studies Alexa nder] 
et al.| |2QQ8[), the averaged SFR for SGMs would stilllDe 
higher than that of 24 /im sources by a factor of (1.5-2). 
The inferred values for SFE and gas depletion time scale 
TSF range from 69-388 M-^((S'F£;) = 330 M"^) 
and 17-96 Myr {{tsf) = 38Myrs) respectively. For com- 
parison, star formation efficiency (SF E) is typically 



(180 zb 16O)L0Mq^ for local ULRIG s ([Solomon et al 
1997|) and 450^170 L qM~^ for SMGs ("Grev e et al.|2QQ5 



Tacconi et al.||2006p. And the gas depletion time scale 
Tsf ranges from 10-lOOMyrs for SMG , and a factor of 



10 longer for loc al LIRGs and ULIRGs (Solomon & Van- 



den Bout 



2005). Figure [5] displays the Lfir/L^q ratio 



versus L^q (or equivalently cold molecular gas), visu- 
ally illustrating our conclusion that these bright 24 /im 
selected z ^2 ULIRGs have less cold molecular gas and 
smaller star formation efficiency than that of SMGs. 

3.3. CO velocity widths and dynamical mass estimates 

The CO velocity width (FWHM) can be measured by 
using either a single or double gaussian fit or directly 
from the data (full width at the half peak intensity). 
Table [4] shows the results using both methods. Assum- 
ing a single component and without profile fitting (Ta- 
ble [4] FWHM2), we have velocity widths ranging from 
128Km/s to 824km/s, with a median value of 406km/s. 
If we use two velocity components for sources with dou- 
ble peak profiles, the averaged velocity width reduces 
to 275 km/s. With a single component fitting. Figure [6] 
compares the CO velocity widths of our sources with that 
of SMGs and QSOs. Here the SMG comparison samp le is 
listed in Table [6] for 14 objects compiled from , Greve et al.| 
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TABLE 6 
Comparison Sample I SMG data 



ID 






T " 

CO(3 — 2) 


^ FWHM 


double 


Ref. 








K km/s pc^ 


km/s 






SMMJ023956-0134 


1.0620 


8.36e+12 


2.1e+10 


780 ± 60 


yes 


(1) 


SMMJ023951-0136 


2.8076 


1.82e+13 


4.8e+10 


1360 ± 50 


yes 


(1) 


SMMJ044307+0210 


2.5094 


3.04e+12 


l.le+10 


350 ± 60 


yes 


(1) 


SMMJ094303+4700 


3.3460 


1.67e+13 


2.7e+10 


420 ± 50 


yes 


(1) 


SMMJ123549+6215(HDF76) 


2.2021 


8.90e+12 


5.0e+10 


600 ± 50 


yes 


(2) 


SMMJ123707+6214(HDF242) 


2.490 


4.00e+12 


2.0e+10 


430 ± 60 


yes 


(2) 


SMMJ131201+4242 


3.408 


1.18e+13 


5.3e+10 


530 ± 50 


no 


(1) 


SMMJ140103+0252 


2.5653 


5.51e+12 


1.7e+10 


200 ± 40 


no 


(1) 


SMMJ163554+6612 


2.5168 


1.33e+12 


3.7e+9 


500 ± 100 


yes 


(1) 


SMMJ163650+4057(N2850.4) 


2.3853 


1.56e+13 


7.0e+10 


840 ± 110 


yes 


(2) 


SMMJ163658+4105(N2850.2) 


2.4500 


2.03e+13 


5.7e+10 


870 ± 80 


yes 


(2) 


SMMJ163706+4053 


2.3800 


2.00e+13 


3.0e+10 


830 ± 130 


no 


(1) 


EROJ164502+4626 


1.4950 


9.69e+12 


4.1e+10 


400 ± 20 


no 


(1) 


SMMJ221735+0015 


3.0990 


1.20e+13 


3.8e+10 


780 ± 100 


no 


(1) 



: 1 and 



^ All [Greve et al.| (|2005|); [Tacconi et al.| (|2006| |2008|) assume that Tco(3-2)/rco(i- 
a = O.8M0/[K,km/spc'^] for the molecular gas mass calculation. 
^ The line velocity width is from a single gaussian fit to the spectral data. 

^ References: (1) Greve et al. 2005, MNRAS, 359, 1165. This paper gives the original references 
for some of the sources which were not observed by Greve et al. 2005. (2) Tacconi et al. 2008, 
ApJ, 680, 246 
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Fig. 5. The plot 

L^q[K km/s pc^] versus L^^ 

infrared and ultra-luminous galaxies from [Solomon et al.| (fl997) 
the red, open triangles are our data points, t lie blue, open squares 
the SMGs the published data from Greve et al. (2005); Tacconi et 
[al. (2006. ^2008 j , and the green crosses the ^(jSO <fe radio galaxy 
oala compiled by [Solomon ' 



Vanden Bo^ ^005\ 



( |2QQ5l ) ; [Tacconi et al | f2QQ6, 2008). The QSO comparison 
sample is compiled by Solomon fc Vanden Bout (2005) 



(also see Greve et al. 



( 2005| j. The right panel presents 
the velocity distributions, with dashed lines marking the 
median values of 406, 635 and 350 km/s for Spitzer z^2 
ULIRGs, SMGs and QSOs respectively. As the local ref- 
erence, the z ^0 U LIRGs have a median CO velocity 



width of 240 km/s ( [Sanders et al.||1991l 
1997[). Kolmogorov-Smirnov tests on ti 



Solomon et al. 
ese three small 
samples produce the probabilities of 0.23 and 0.84 respec- 
tively, for AV(FWHM) of 24 /im ULIRGs versus SMGs 
and 24 /im ULIRGs versus QSOs drawn from the same 



parent distributions. 

Recent N-body/Smoothed-Particle Hydrodynamic 
(SPH) numerical simulation c ombined with 3D poly- 
chromatic radiative transfer by [Narayanan et al. (2009) 
show that SMGs are mostly produced by major mergers 
of gas rich, star forming dominated progenitors, whereas 
24 /im bright ULIRGs at z ~ 2 could be formed by both 
major and minor mergers. This model predicts that CO 
velocity widths of bright 24 jam ULIRGs spread over 
(100-1500) km/s, overlapping with that observed among 
SMGs, but the median value is smaller than that of 
SMGs. 

Spatial distribution of molecular gas holds the 
imprints of kinematic state of a galaxy. Al- 
though the majority of our sample is unresolved, 
two targets, MIPS506 and MIPS16144, show spa- 
tially extended CO emission across regions of ~ 100 
and 46kpc, significant after taking into account of 
their beam sizes of 5.02 x4.37 ^41kpcx36kpc and 
2.44"x2.15" -20kpcxl8kpc respectively. MIPS506 
shows two distinct components with a separation of 
45kpc. The CO spectra in Figure [3] show that MIPS506, 
MIPS16144, and MIPS8342 have double peaked veloc- 
ity profiles, with single and double gaussian fit (solid 
and dashed lines respectively) to the observed data (his- 
togram). A double peak CO velocity spectrum is in- 
dicative of either a rotating disks or of two galaxies in 
an ongoing merger. Our CO data measures 3 out of 9 
(33±19%) sources having double velocity components. 
The low resolution CO data for the SMG comparison 
sample in Table[6]shows 64±21% (9/14) having double or 
multiple velocity components. If we consider pure Pois- 
son statistical errors, these two fractions are consistent 
within la. 

One important question is whether the two velocity 
components in each of the three double peak sources 
are spatially separated. Figure [7| shows the CO con- 
tours overlaid on the 24 /im (pink) and HST H-band 
(black) images, with solid and dashed contours indicat- 
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Fig. 6. — The plot compares the CO velocity width (FWHM) ver- 
sus molecular gas masses of our sources with that of high-z QSO &; 
radio galaxies and SMGs. All velocity widths are based on a single 
cornponent fitting. The SMG comparison sample is listed in Ta- 
ble [m compiled from Greve et al. (2005); Tacconi et al. (2006), and 
the <QSO &; radio galaxy data are trom the compilation ot [Solomon 
fc; Vanden Bout (2005); Greve et al. (2005). The right panel shows 
The velocity distributions ot three populations with the dashed lines 
marking the median values. The velocity width distributions show 
tentative evidence that 24 /im ULIRGs are similar to QSOs, having 
smaller CO velocity width than the averaged value for SMGs. 



ing two separate velocity peaks. Here, the integrated 
CO maps are summed over (4-5) velocity channels (see 
Table [2| for the channel width), covering the respective 
velocity peak. For example, the integrated CO maps 
for MIPS16144 shown in Figure [7| are summed from - 
730km/s to -330km/s for the velocity peak at -605km/s 
and from -330km/s to -70 /km/s for the velocity peak at 
-138km/s. The similar maps are made for MIPS506 and 
MIPS8342, and shown in Figure[7j The rms a values used 
for the contours are in Table |2] The two velocity peaks in 
MIPS 16144 have a clear spatial separation of 1.3 , which 
corresponds to 10.9 kpc. The two velocity components in 
MIPS8342 are spatially too close to be resolved. In the 
case of MIPS506, the separation is ~ 5.5 , correspond- 
ing to 45 kpc. Such a large distance, particularly for cold 
molecular gas, suggests that MIPS506 is probably a pair 
of merging galaxies, with substantial dust extinction in 
the rest-frame optical band. This explains the very faint 
multiple blobs (2.8 ~ 24 kpc separation) in the H-band 
image (6 x 5 ) shown in Figure [sj We have examined 
our shallow 3.5- 8 /im IRAC images to see if the galaxy 
pair detected in the CO map is also detected in other 
bands, only an unresolved object is detected at 3.5 and 
4.8 /xm, and no detections at other two IRAC bands. The 
physical model for MIPS16144 and MIPS8342 is prob- 
ably also merger, although rotating disk model is also 
consistent. Additional supporting evidence for Spitzer 
ULIRGs being mergers comes from the HST/NICMOS 
H-band images, with the detailed discussion in §3.4[ 

For MIPS16144 and MIPS506 with spatially re- 
solved CO maps, assuming a two-body interaction, 
we can approximate the total dynamical mass with 
Mdyn = R*v^irc/G = R*AVV(sin2(i)G), where sin(i) is 



TABLE 7 
Comparison Sample II QSO data"^ 



ID 


^co 


Mbh 
109 

Mo 


Mgas 

IQio 

Mo 


Mdyn 

10^0 

Mo 


incl. angle 


APM08279+5255 


3.91 


23 


13 


22 


25deg 


PSSJ2322+1944 


4.12 


1.5 


1.7 


4.4 


no corr. 


BRI1335-0417 


4.41 


6 


9.2 


10 


no corr. 


SDSSJ1148+5251 


6.42 


3 


2.2 


5.5 


65deg 



^ The black hole masses are derived from UV spectroscopy, the 
gas and dynamic masses from the high spatial resolution CO 
observations. References: [Walter et al.| ( |2004| ; [Riechers et al.| 
( |2008ab, 2009; . 



the inclination angle, R is the separation between two 
merging objects, and AV is the observed CO peak ve- 
locity difference. R and AV are (10.9 kpc, 463 km/s) 
and (45kpc,174km/s) for MIPS16144 and MIPS506 re- 
spectively. We derived Mdyn = 5.4x 10^^sin~^(i) and 
3.2xlO^^sin-2(i)M0 for MIPS16144 and MIPS506. If 
assuming an averaged value of (sin^(i)) = 2/3, the gas 
fraction, Mgas/Mdyn, is 4% and 3% for MIPS16144 and 
MIPS506 respectively. Although several recent papers 
have published gas fractions for various types of galax- 
ies, we caution that the equation used for calculating 
the dynamical masses based on observed e mission line 
FWHM could differ by a factor of (2-5 ) 



see 



(I2QQ3) 



et al. 



Erb et al. (2006); Tacconi et al.|(2006); Riechers 



Neri et al. 



TEe stellar half-light radii have been measured from 
the HST/NICMOS H-band (1.6 /im) images for our sam- 
ple dPasyra et al ||2008t |Zamojski et al.||2009p . Us- 
ing Mdyn = Ri/2*vL/G = Ri/2nAV--HM)7pn2(i)G) 
and assuming that cold molecular gas disks are smaller 
than stellar half-light radii, we can also set limits on the 
dynamical masses of additional 5 sources in our sam- 
ple. This assumption is usually correct f or nearby galax- 
ies dSanders et aTjllOQll [Solomon et al.| |1997). We note 
that our sources probably have fairly high dust extinction 
at the observed H-band (rest-frame 5000A), which will 
cause systematic under-estimates of the intrinsic stellar 
sizes. Using the half-light radii are listed in Table [5) 
we obtained Mdyn (0.2-2.2) x 10^^ for MIPS8327, 
MIPS8432, MIPS15949, MIPS16059 and MIPS16080, us- 
ing an averaged inclination angle. These rough estimates 
are on the same order of magnitude in comparison with 
the averaged dynamical mass o f SMGs, 1.2xlO^^M0 
dOreve et al]|2005[ [Tacconi et ar||2Q06^ . 



3.4. Our targets are likely mergers 

The CO images and spectra in Figure [T] and Fig- 
ure [3] have revealed that at least one source, MIPS506, 
is a merger involving two galaxies with substantial 
cold molecular gas. In addition, MIPS8327 and 
MIPS 16059 are likely interacting galaxies because their 
HST/NICMOS H-band images show close companions 

within sepa rations of 3.3 a n d 4 k pc (0.4 & 0.5 ) re- 
spectively. "Zamojski et alj (2009) performed detailed 
2-dimensional light profile fittmg to the HST images, as- 
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same ratio for these ^^3 Q SOs is on the order of 1:10 



Fig. 7. Overlaid on the 24 /im (pink) and NICMOS H-band 
(black) images are the CO contours of the integrated maps over 
the positive (solid lines) and negative (dashed lines) velocity peaks 
for the 3 sources with double peak profiles. The solid and dashed 
contours are for the CO velocity peaks at (-138, -605)km/s, 
(160,-15) km/s and (-250,-520) km/s for MIPS16144,MIPS506 and 
MIPS8342 respectively. The contour starts at la w ith a step of la 
for all three panels. For the detail, see the text in §3.3| 



suming a sersic profile modeP] Figure [s] compares the 
original H-band images with the model subtracted, resid- 
ual images (reversed intensity scale) for 5 of our sources 
with CO detections, MIPS506, MIPS16144, MIPS8342, 
MIPS16080 and MIPS15949. Multiple blobs an clear 
tidal tails in Figure|8]suggests that MIPS506, MIPS16144 
and MIPS 16080 are probably the products of galaxy 
interactions. Although the original H-band images of 
MIPS8342 and MIPS15949 show isolated early type 
galaxies, but the faint disks/rings in the residual images 
indicate that they could be late stage mergers. Majority 
sources in our CO sample have merger signatures either 
in their CO maps or in the H-band images. Consider- 
ing the published studies of SMGs based on HST im- 
ages and integral spectroscopy ( Chapman ^t al. 2003; 
Conselice et al. 2003 Swinbank 'et al. 2006|, we spec- 
ulate that SMGs, bright 24 /xm ULIRGs, and gas rich 
QSOs are connected by an evolutionary model involving 
mergers, with SMGs being in an early stage of merg- 
ers of two gas rich, star forming progenitors, and bright 
24 /im ULIRGs with substantial obscured AGN in a later 
merger stage, while QSOs are fully merged with massive 
black holes in place, producing powerful feedback clear- 
ing away dust obscuration. This simple picture does not 
explain how SMGs and gas rich QSOs could have such a 
different history of star formation and black hole growth. 
Specifically, most SMGs have very weak AGNs, and their 
estimated black hol e to stellar mass ratio is roughlv on 



the o rder of 1:10^ (Borys et al. 2005 



2008|, whereas, recent studies (see 



5||^ 

TaK 



Ale2cander et al 
le [7| found the" 



4 /(r) = Ie*exp(— * [(r/re)-*^/^ — 1]), is a function of sersic 
index n. Sersic index n is 4 for de Vaucouleurs profile and 1 for 
the exponential profile. 



( Riechers et al.|200'8a , 2009). For our objects, black hole 
mass is calculated from the AGN bolometric luminosity 
estimated in the SED decomposition (see Figure [2]). Stel- 
lar mass is computed from the rest-frame Lh assuming 
a Mstars/Ln ^ 0.48 from a 2 Gyr SSP template of Maras- 
ton (2005). The inferred Mbh/M 

stars is on the order of 
1:10^. Again, this suggests that Spitzer ULIRGs may be 
a transitional type of sources between SMGs and QSOs. 
To really understand the physical relation among these 
three types of sources, we will need better UV optical 
spectra, complete SEDs, and high resolution CO data 
for Spitzer ULIRGs. 

4. SUMMARY AND DISCUSSIONS 

We report CO interferometric observations of nine 
z^2 ULIRGs with 

S24fim ^ 1 mJy. Combining high spa- 
tial resolution HST images, Spitzer mid-IR spectra, pho- 
tometry and the CO data, we find the following results: 

• Of the nine sources observed at PdBI, eight 
sources have significant detections of the CO J(3- 
2) or J(2-l) transition with intensity Ic ~(0.5- 
1.5)Jykms~^ at (4.7-9.7)(T. CO line emission 
is detected from sources with a variety of mid-IR 
spectral type, including strong PAH, deep silicate 
absorptions and mid-IR power-law. The observed 
CO J(3-2) [or (2-1)] line luminosities L^q are (1.3- 
5.5)xl0^^ Kkm/spc^, yielding cold molecular gas 
masses of (1 -4.4) x IO^^Mq, based on the assump- 
tion of Tco(3-2)/Tco(i-o) = l and the CO (1-0) 
luminosity to H2 gas mass conversion factor of 
0.8 M0/[K Km/s pc^]. Based on small statistical 
samples, we find tentative evidence suggesting that 
Spitzer ULIRGs have a factor of 2 less cold molec- 
ular gas than what observed among SMGs and gas 
rich QSOs. 

• The CO velocity width distribution of Spitzer 
ULIRGs is similar to that of QSOs with CO de- 
tections, but on average a factor of 1.5 smaller 
than that of SMGs, with caveat that the avail- 
able CO statistics is small. Two objects, MIPS506 
and MIPS16144, have spatially resolved CO emis- 
sion across regions of 100 and 46kpc respectively, 
significantly extended after taking into account of 
their beam sizes and shapes. Three of our nine tar- 
gets (33%) have double peaked CO spectra, in com- 
parison with 64% observed among SMGs. Two of 
these three, MIPS506 and MIPS16144, have their 
two velocity peaks spatially separated by ~ 45 kpc 
and 10.9 kpc respectively. Such a large spatial sep- 
aration of two distinct CO knots suggests that 
MIPS506 is a pair of merging galaxies. The in- 
ferred dynamical masses are 5.4x 10^^ sin~^(i) and 
3.2xlOiisin-2(i)M0 for MIPS16144 and MIPS506 
respectively, yielding the gas fraction (Mgas/Mdyn) 
of 8% and 6% if using an averaged inclination an- 
gle. 

• Together with spatially resolved CO data, our 
HST/NICMOS H-band images suggest that ma- 
jority of our sample have signatures of mergers at 
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Fig. 8. — The HST H-band images and their residual images with , 
distinct galaxy pairs. The intensity scale is reversed for the residual 
morphological signatures suggesting dynamical interactions. 

various stages of dynamical interactions, includ- 
ing close companions, tidal tails, and rings in the 
residual images. We hypothesis that SMGs, bright 
24 jam ULIRGs and gas rich QSOs may be con- 
nected by the same evolutionary model in which 
the formation mechanism is primarily mergers, 
with SMGs mostly at an early passage of two gas 
rich, star forming galaxies, bright 24 /am ULIRGs at 
a later stage of dynamical interaction, with gas and 
dust surrounding the central black holes, resulting 
smaller observed CO velocity width, and finally, 
the gas rich QSOs representing the final merged 
stage. 



a sersic profile model subtracted for the 5 sources which do not show 
images in order to bring out the faint features. These images reveal 
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